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Abstract 

The field-induced quantum criticality of compounds with ferromagnetically coupled 
structural spin units (as dimers and ladders) is explored by applying Wilson's renor- 
malization group framework to an appropriate effective action. We determine the 
low-temperature phase boundary and the behavior of relevant quantities decreasing 
the temperature with the applied magnetic field fixed at its quantum critical point 
value. In this context, a plausible interpretation of some recent experimental results 
is also suggested. 
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The study of quantum phase transitions (QPT's) and low-temperature prop- 
erties close to a quantum critical point (QCP) of a wide variety of materials 
has recently attracted considerable attention and constitutes today a topical 
subject in condensed matter physics [1]. 

The main methods of tuning a system toward a QCP are essentially based 
on manipulation of doping, pressure and magnetic field. A lot of experiments 
on spin compounds has shown that the magnetic field is the most convenient 
non-thermal parameter to control the distance from a QCP. 

Emerging magnetic field-induced QPT's have been observed in quantum anti- 
ferromagnetic (AFM) compounds as KCuCh [2], TlCuCk [2,3], BaCuSi20e 
[4], and Cu2{Cr,HY2.N2)2ClA [5] which consist of weakly coupled low-dimensional 
structural spin units such as dimers [2,3,4] or ladders [5,6]. In absence of a mag- 
netic field, these systems exhibit a gap between a singlet ground state and the 
lowest triplet excitation. When the field is switched on and increased, the gap 
decreases by Zeeman effect and a field-induced QPT occurs at a critical field 
which measures the amplitude of the original gap. 

Recently, also the case of weakly ferromagnetically coupled spin units, al- 
though less explored then the AFM one, has attracted a lot of interest [7,8,9,10,11,12] 
motivated by the possible existence of dimer and ladder materials with effec- 
tive ferromagnetic (FM) interactions between the basic spin units [7,9,12,13,14]. 
Inelastic neutron scattering (INS) investigations on the spin dimer compounds 
Cs2,Cr2Br^ [7], on the layered cuprate SruCu2404i [12] and on CaV205 [11], 
which contains layers of coupled two-leg ladders, appear particularly meaning- 
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ful in this direction. An intriguing feature is that these materials are charac- 
terized by frustrated inter-units microscopic coupUngs. Nevertheless, the INS 
predictions suggest that their essential physics can be captured by means of a 
spin model where the original frustrated inter-units couplings are replaced by 
FM ones. A zero-temperature numerical study of this effective simplified spin 
model has been performed in Ref. [9] and the results are in agreement with 
previous conventional investigations [15]. 

Effective intcr-units FM couplings may also arise from different mechanisms. 
An important example has been considered in Ref. [11] where a unified picture 
of recent INS susceptibility data [16] for stripe-ordered Lai^/gBai/sCuO^ is 
given based on a model of two-legs spin ladders where an effective FM inter- 
ladder coupling results from integrating out the degrees of freedom in the 
stripes. 

Support to the ferromagnetically coupled structural spin units scenario is given 
also by first-principles and Monte Carlo calculations of electronic structures 
about the nature of the spin-singlet ground state in the dimer compound 
CaCuGe20e [14]. The results obtained for the susceptibility and magnetiza- 
tion behavior give indeed evidence of effective FM inter-dimer couplings and 
of long-range interactions effects. In addition there are indications [9] of com- 
pounds which consist of ferromagnetically coupled ladders, as SrCu20^ [17]. 

To gain insight into the properties of spin materials with weak FM inter- 
units couplings, also in the presence of a magnetic field, the original com- 
plex spin Hamiltonian is conveniently mapped in an effective spin- 1/2 FM 
XXZ model with spins, exchange couplings and longitudinal magnetic field 
expressed perturbatively as linear combinations of the original microscopic 
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ones [18,7,8,9,10]. The final Hamiltonian has the general structure 

N N 

h^-Y: [MS!s^ + sfs]) + K,jS^s^ -hY^st (1) 

i,j=l 1=1 

where S°'{a = x,y,z) denote the effective spin components at site i (which 
is a spin unit index in the original spin model) of a lattice with N sites and 
Jij > 0, Kij > (with Jii = Kii = 0). 

A systematic study of the thermodynamics of weakly ferromagnetically cou- 
pled spin units was performed [18] more than three decades ago within a 
mean field approximation (MFA). Further few studies for specific problems at 
zero temperature and, in some cases, in absence of a magnetic field, have been 
achieved at numerical, Hartree-Fock and MFA levels [7,8,9,10,14,18]. However, 
reliable experimental and theoretical studies of the magnetic properties and of 
the phase boundary of these systems close to the field-induced QCP are still 
lacking at the present time. 

The aim of this letter is to give a contribution in this direction by using 
an appropriate functional representation of the FM spin model (1) [19,20] 
and a wilsonian renormalization group (RG) approach, which is the most 
valid and reliable tool to take properly into account fluctuations effects close 
to quantum and classical critical points. For generality purposes and with 
the intent to interpret some experimental findings [21,22], we refer to a d- 
dimensional spin lattice and include the possibihty of long-range FM couphngs 
which decrease with the distance between the spins as a power law of the type 
with (7 < 2. The value a — 2 corresponds to nearest- neighbor FM 
spin-spin interactions. 

The action appropriate to describe the iow-tempcraturc properties of the FM 
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model (1) was derived almost three decades ago [19] and sounds as 



^ Y.U{uo,Vo;{q.})r{qi)r{q2mq3)^{q4) (2) 



where g = (^,0;^ a,ndU {uo,Vo; {q^}) = %^+^^.if3+^^K(^a;,j+a;,2;c.i3+a.U+f^0(^a;,j;a;,3(^a;,,;a;,J. 

Here, ip{k,uji) is a complex field related to the in-plane magnetization, k de- 
note the wave vectors (we assume a cut-off A = 1 related to the original spin 
lattice), oji — 2t^ITq (I — 0, ±1, ±2, ...) are the bosonic Matsubara frequencies 
and V is the volume of the system. In Eq. (2) Tq (xT {T is the physical tem- 
perature), ro (X h — [J(0) — K{0)], vo oc K{0), where J(0) and K{0) are the 

— * 

{k — 0)-Fourier transforms of Jij and Kij in Eq. (1). The explicit expressions 
of To and ro, uq, vq in terms of the coupling parameters in the Hamiltonian (1), 
inessential for present purposes, can be found in Ref . [19] . Of course, if we put 
Vq = {Kij = in Eq. (1) ), d — 3 and a — 2, the action (2) reduces to the 
well-known one for a 3-dimensional XY model in a transverse field [23] and 
for a dilute gas of hard-core bosons with chemical potential // = — ro, which is 
believed [24,25,26] to be appropriate for a description of the dimer compounds 
XCuCh {X = K, Tl) close to their QCP's. 



Applying the Wilson RG transformation to the action (2) and working to 
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one-loop approximation we obtain the equations 



dr , Kd 



{v + 2u) Fi (r, T) 



1 + r. 

^=i^-d)u-^u^ (r, T) + F3 (r, T)] - ^^^.^^ 



(3) 



2(7(l + r)V 

for which the appropriate Fisher exponent is 7; = 2 — cr. In Eqs. (3), / de- 
notes the RG rescahng parameter, Fi (r, T) = (1/2) coth[(l + r) /T] , F2(r, T) = 
-dFi (r, T) /ar and F3 (r, T) = Fi (r, T) / ( 1 + r ) . 

The (To = 0)-critical properties can be simply obtained setting T = in 
Eqs. (3). For case of interest d > a, when the quantum gaussian fixed point is 
stable, one finds MFA criticality in terms of (ro — Tqc) oc {h — he) as Tq ^ rj,, 
where Tqc = —{Kd/4:d){vo + 2uq) localizes the QCP. Of course, the corre- 
sponding (To = 0)-critical value he of h can be easily found from the explicit 
expressions of Tq, Mq, Vq in terms of the original coupling parameters. 

The wilsonian classical critical regime at finite temperature for d < 2a can be 
analyzed defining u — uT, v — vT, and then setting T{1) — > 00 as Z — > 00 in 
the RG equations for the new coupling parameters u and v. 

Here, we are interested to solve Eqs. (3) ior a < d < 2a hmiting ourselves 
to quantum or classical gaussian regime. This can be performed in the low- 
temperature limit and working to leading order in the coupling parameters. 

With u{l) ~ Uoe-^'^-''^\ v{l) ~ Voe-^'^-"^^ and T{1) ~ 6"% {z ^ a is the 
appropriate dynamical critical exponent), r(/) is obtained through the non- 
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linear relevant scaling field g{l) — r{l) + {Kci/M){v{l) + 2u{l)) which scales 
as 



d+a 



+ (^0 + 2t.o) T,^ / dy^-, , (4) 

with go^ro- roc- 

We now stop the renormalization procedure at a scale T » 1 to be determined 
setting g{l*) — 1. Then, from Eq. (4), with Tq -C 1 but arbitrary T(/*) = 
e'^'*To, we have for the dimensionless inverse susceptibility x = (x/xo)~^ = 
e~°''* , the self-consistent equation 



+ ^(y, + 2u„)lt" r""dy^. (5) 
Za Jo Qy — \ 



which contains all the relevant physical information and allows us to explore 
the full quantum critical region in the phase diagram for d > cr. 

The phase boundary equation roc(To) in the (ro, To)-plane, ending in the QCP, 
can be obtained setting x = (Z* = oo) in Eq. (5). One obtains 



oc 



goc {Tq) = Toe (To) - ro 



vo + 2uo) Tt, (6) 



2{d-a) 2a 

with Ad^u = T{d/a)({d/a). It is worth noting that, if one assumes vq = 0, we 
have roc{To)-roc oc l^Hc{T) = hc{T)-hc oc T'^/'" and hence the result i/j = d/a 
for the phase boundary or shift exponent defined as AHc{T) oc T^. For d — 3 
and cr = 2, we find, as expected [24,25], ijj = 3/2. In contrast, with f o 7^ 
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and for sufficiently low temperature, Eq. (6) yields ip — 1 which is a universal 
value independent of d and the interaction parameter a. Thus, the structure of 
the phase boundary equation (6) close to the QCP suggests the existence of a 
crossover temperature Tq between two different regimes characterized by the 
values d/cr and 1 of the exponent 'ip decreasing the temperature. By inspection 
of Eq. (6), one immediately sees that (with small — coupling) 



T* 



avo 



{d - a) Ad^cr {vo + 2uq) 



[X(0)]^, 



(7) 



and the mentioned crossover is described by the effective exponent 



d — a 

T o- 



d\n.T 



1+T- 



(8) 



with T = To/To* and 1 < ^/'e//(r) < d/a. For case d = 3 and a = 2, Eq. (8) 
yields V^e// (r) = [l + (3/2)rV2] / (i + t^I^^ with 1 < V^e// < 3/2. 

Although all the quantum critical properties in the gaussian region close to 
the QCP can be obtained solving the self-consistent Eq. (5), we hmit ourselves 

to the experimentally relevant case tq = tqc {h = he) as Tq — > 0. A solution 
exists for Toe"^* = T{1*) > 1 and we find 



— 2(d-<T)^0 ^ 



2(T 



KdVp 
2{d-a) 



To, To < T* 



2(7 



t;o + 2mo) To-- , To>To*. 



(9) 



This implies that also for susceptibility x ~ e'^'* and correlation length ^ ~ 
T^i/^ ~ e'* one can define effective exponents yielding 7e// (r) = -^e// (t) and 
Ueff (r) = (l/(7)'0e// {'t)- For the realistic values d — 3 and a — 2, one has 

(X ~ To-\c ~ To"'/') and (x ~ To-'/',^ ~ To-'/') for To < T* and To > T*, 
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respectively. The Iow-Tq behaviors of other quantities can be obtained using 
the standard RG machinery. 

The previous results, and in particular the phase boundary exponent ip (or, 
more properly, ipeff{T) in our RG scenario), may constitute a good basis for 
comparisons with available and possible future experimental data. Indeed our 
RG analysis suggests that, close to the QCP, one must expect boson-like be- 
havior (-0 = d/a) only above a characteristic temperature T* oc Tq{vo,Uo) 
depending on the coupling parameters in the original effective FM XXZ spin 
model. Decreasing the physical temperature below T*, a crossover takes place 
to a different regime where the phase boundary is characterized by a different 
universal exponent ip — 1. Oi course, the amplitude of the crossover region 
reduces to zero decreasing vq. In particular, when vq = or negligibly small, 
only the exponent ip — d/a should be observable. In contrast, if the crossover 
region is sufficiently large, only the linear behaviour in T is expected at suffi- 
ciently low temperature. 

A situation of this type seems to occur for the doped dimer compound Tli^^KxCuClz 
for which the effect of randomness on field induced magnetic ordering has been 
recently investigated through low-temperature specific heat measurements for 
potassium concentration in the interval < x < 0.22 [27]. The experimental 
predictions for the phase boundary exponent ■0, relevant for our purposes, can 
be summarized as follows. 

First, a ree valuation of -0 for a; = yields the value ■0 = 1.67 which is close to 
V'bec = 3/2 derived, for compound TlCuCls, from the theory of Bose- Einstein 
condensation of triplet excitations. 

Next, a systematic study for a; 7^ shows that the disorder produces a sensible 
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change in the critical scenario. In particular, the phase boundary is accurately 
measured for temperature below 2 K and the relevant features are: 

i) the exponent decreases systematically increasing x from the {x = 0)-value 
■0 = 1.67 and tends to ■0 = 1 for x > 0.1. This corresponds to an enlargement 
of the phase boundaries obtained for different values of x; 

ii) the phase boundary observed for x > 0.1 is almost a linear function of 
temperature. These findings are clearly shown in Fig. 4 of Ref . [27] . 

To explain these results, it was conjectured in Ref. [27] that a Bose glass of 
triplons, in the sense argued by Fisher et al. [28], could be an appropriate 
model to describe properly the quantum criticality of Tli^^KxCuCl^ with 
x 7^ 0. Nevertheless, the Bose glass theory, as applied to this compound close 
to the QCP, is expected to predict [27] a phase boundary exponent ■0 < 1/2 as 
T — > 0, in drastic disagreement with the available measurements. To solve this 
puzzle it was also suggested that a crossover from the convex form (■^ > 1) to 
the concave form (V' < 1) of the phase boundary should occur very close to 
the QCP. However, this type of crossover is not observed in the experiment 
and hence the problem remains an open question. 

The results here obtained allow us to speculate that the enlargement of the 
phase boundaries and the almost linear dependence of the critical field on 
temperature for x > 0.1 [27] can be simply interpreted in terms of the previous 
RG scenario involving the characteristic crossover temperature T* ~ Vq^^'^ '^^ 
with a — 2 and d — 3. To show this, we adopt for Tli^^KxCuCl^ the idea 
that the doping, as frustration, induces an effective FM coupling between spin 
dimers. Furthermore we assume, reasonably, that vq, and hence T*, vanishes or 
its effect becomes negligible when x decreases to zero (this hes on the feature 
that the action (2), with vq = 0, appears to work for the dimer compound 
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TlCuCh [24,25]). 

For X — Q our analysis predicts that the phase boundary behavior for T — > 
is characterized by the Bose exponent = 3/2, as expected for TlCuCl^ [26]. 
For X ^ 0, the effective FM — coupling in the Hamiltonian (1) becomes 
active and a crossover is expected to occur from ■0 = 3/2 to ■^^l decreas- 
ing the temperature to zero across T*. When a; < 0.1 [27], the "temperature 
window" (0,7*) starts to be experimentally accessible and a -0 < 3/2 is mea- 
sured. Increasing x (above 0.1), this window becomes larger and larger, the 
temperature can be decreased to zero working below T* and one recovers the 
exponent ijj — 1 which characterizes the linear temperature behavior of the 
phase boundary observed in the experiment close to the QCP. Besides, since 
for T ^ it is found 3/2 > > 1 increasing x, one can argue that an en- 
largement of the phase boundaries takes place, again in agreement with the 
experimental data [27]. 

From our RG predictions one can obtain also some insight into the x-dcpcndcnce 
of the crossover temperature T* as given by Eq. (7) in terms of the effective 
longitudinal spin couphng K{0), with vq oc K{0) [19]. Indeed, since one expects 
that = for a; = [24,25], we can assume plausibly that ^0(2;) oc K{0) ~ x 
for small potassium concentration. Hence, Eq. (7) yields T* ~ and, for 

d — 3 and a — 2, one has T*(x) ~ x^, result which could be tested experimen- 
tally. From a theoretical point of view, these results should be extracted from 
an appropriately chosen effective spin dimer model and then using the explicit 
relations between the bare couphngs and the effective FM ones obtained by 
mapping [18] the original Hamiltonian in the FM XXZ model (1). A study 
of this not easy problem is in progress and we plan to present a quantitative 
scenario in a future work. 
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In any case, further reliable experimental and theoretical studies on ferro- 
magnetically coupled structural spin units are desirable and we hope that this 
preliminary contribution will stimulate more intensive research activity on the 
subject. 
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